Canonical quantization, path integral representations, and pseudoclassical description 
of massive Weyl neutrinos in external backgrounds 
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We study massive 1/2-spin particles in various external backgrounds keeping in mind applications 
to neutrino physics. We are mainly interested in massive Majorana (Weyl) fields. However, massive 
neutral Dirac particles are also considered. We formulate classical Lagrangian theory of the massive 
Weyl field in terms of Grassmann-odd two-component spinors. Then we construct the Hamiltonian 
formulation of such a theory, which turns out to be a theory with second-class constraints. Using 
this formulation we canonically quantize the massive free Weyl field. We derive propagators of the 
Weyl field and relate them to the propagator of a massive Dirac particle. We also study the massive 
Weyl particles propagating in the background mater. We find the path integral representation for 
the propagator of such a field, as well as the corresponding pseudoclassical particle action. The 
massless limit of the Weyl field interacting with the matter is considered and compared with results 
of other works. Finally, the path integral representation for the propagator of the neutral massive 
Dirac particle with an anomalous magnetic moment moving in the background matter and external 
electromagnetic field, as well as the corresponding pseudoclassical particle action are constructed. 
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I. INTRODUCTION 

The study of neutrinos is one of the fastest developing 
areas of elementary particle physics. This is mainly due 
to recent experimental achievements in measuring the pa- 
rameters of the neutrino mixing matrix [l[ . Great experi- 
mental efforts have been devoted to clarify the question of 
whether neutrinos are Dirac or Majorana particles (see, 
e.g., Ref. 0). Nevertheless the problem of the neutrino's 
nature still remains open. It should be noted that the 
Majorana type of neutrino masses are predicted in some 
models of the neutrino mass generation, like the see-saw 
mechanism Q. However, one cannot exclude the possibil- 
ity that neutrinos will turn out to be Dirac particles. In 
the present work, we study, in the main (except Sec. |V| . 
Majorana (or Weyl) fields representing neutrinos. 

The first quantum field theory treatment of massive 
Majorana neutrinos was made in Ref. (f|, where it was 
proposed to describe a massive Majorana field in terms 
of two-component Weyl spinors. Some processes in- 
volving neutrinos, like j3 decay, were also discussed in 
Ref. [J] using Weyl field formalism. Mixed massive Ma- 
jorana neutrinos possessing anomalous transition mag- 
netic moments in an external magnetic field were stud- 
ied in Ref. Q. The treatment of massive Majorana 
fermions based of spinors, which belong to the nonstan- 
dard Wigner classes and violate the Lorentz invariance, 
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was proposed in Ref. Q. Recently a detailed quantiza- 
tion of massive Majorana particles, described in terms 
of c-number Weyl fields, in vacuum and in background 
matter, was considered in Ref. Q. 



In the present work, we continue the rigorous study of 
the massive Weyl field in vacuum as well as in background 
matter in the framework of classical and quantum theory. 
In Sec. HH we start with the construction of the classical 
Lagrangian and Hamiltonian formulations of the mas- 
sive Weyl field in a vacuum using the Grassmann- valued 
two-component spinors. Then, based on the constructed 
Hamiltonian formulation, we canonically quantize this 
theory. We also introduce propagators of the massive 
Weyl field and relate them with the corresponding prop- 
agators of the massive Dirac field. The interaction of the 
system of spin- 1/2 fcrmionic fields with general external 
fields is described in Sec. Mil Then, in Sec. lIVl we discuss 
the propagation of a massive Weyl field in a background 
matter. We find the path integral representations for all 
the propagators and derive the corresponding pseudo- 
classical actions; in particular, the pseudoclassical action 
for massive Weyl particles moving in a background mat- 
ter. The massless limit of a Weyl field interacting with 
background matter is also considered. Finally, in Sec. [V] 
we briefly discuss the massive neutral Dirac particles pos- 
sessing an anomalous magnetic moment, moving in back- 
ground matter and in an external electromagnetic field. 
We summarize the obtained results in Sec. IVII 
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II. CLASSICAL AND QUANTUM THEORY OF 
A FREE MASSIVE WEYL FIELD 

As we mentioned in Sec. U the most prominent can- 
didates to be described in terms of Majorana fields are 
neutrinos. Now, it is a well-established fact that there 
are three active neutrino generations (see, e.g., Ref. Q) 
and one has a nonzero mixing between them. In some 
theoretical models a number of sterile neutrinos is also 
predicted @. After the diagonalization of the neutrino 
mass matrix, one gets a number of fermionic fields which 
are Majorana ■ Therefore, if we deal with Majorana 
particles in vacuum, we can study a single massive held 
without loss of generality. Our results can be straight- 
forwardly generalized to include several neutrino genera- 
tions. 

A fermionic spin- 1/2 held ip (x) of the mass m obeys 
the Dirac equation: 

(i7 M 9 M - m)ip (x) = 0, x = (x») = (x°,x) . (2.1) 

For our purposes, it convenient to chose the Dirac 7 ma- 
trices in the chiral representation: 



7 M = 



-a" 



-cr" 

<t» = {o d ,-<t), *» = (a ,*), 



(2.2) 



where a° = <7q is the unit 2x2 matrix, and cr are the 
Pauli matrices (see, e.g., Ref. (Tl|)- 

Let us represent the bispinor ip in terms of two- 
component Weyl spinor r\ as follows: 



102 V 

V 



(2.3) 



Note that such a bispinor satisfies the Majorana condi- 
tion i\) c = ij 2 ip* = ijj. The Dirac equation (12.11) implies 
the following equation for the spinor 77, 



cr^dyj] + ma 2 rj* = 0; 



(2.4) 



see, e.g. Ref. [12j. 

Note that Eq. (|2.4[) is the Euler-Lagrange equation for 
the Lagrangian 



C = vr^a^dyXt - \-m (77 1 a 2 r\ - rfoitf) 



(2.5) 



In this Lagrangian, we treat 77 as Grassmann-odd 
fields [13j . In particular, due to this reason the mass 
term on the right-hand side of Eq. (|2.5j) is not zero. 

In case we study massless spin-1/2 particles, the cor- 
responding bispinor ipo satisfies the equation 



(2.6) 



i^d^o = 0. 



It should be noted that there are solutions of Eq. (|2.6j) 
which satisfy the conditions PriPo = (left-handed 
bispinors) or P^ipo = (right-handed bispinors), where 



Pi 



R,L 



(l± 7 5 )/2, 7 5 =i 7 °7 2 7 2 7 3 



(2.7) 



Considering, e.g., a left-handed bispinor, we can present 
it as 







(2.8) 



Combining Eqs. (|2.6[) and (|2.8|l . we get the wave equation 
for the spinor 770, 



a^dum = 0. 



(2.9) 



The latter equation is the Euler-Lagrange equation for 
the following Lagrangian: 



Co = 'vqla^d^Q. 



(2.10) 



Starting with the Lagrangian in Eq. (|2.5p , we find the 
canonical momenta as 



7T = -^-7- = 1?7 " , 7T 



dfj 



dff 



0, (2.11) 



where the subscript r denotes the right derivatives: see 
Ref. [HI • One can easily see that the system in question 
has two second-class constraints: 

$1 = 7T - irf = 0, $ 2 = 7T*=0. (2.12) 

The extended Hamiltonian reads 

Ui = 7^ + $iAi + A 2 $ 2 , (2.13) 

where 

H = i?7 t (crV)77 + X -m (r) T a 2 r) ~ rfa 2 rf) , (2.14) 

and Ai_2 are Lagrange multipliers which are Grassmann- 
odd fields. These multipliers can be found explicitly, 

Ai = (erV)?? - ma 2 r]*, \ 2 = - (Vrf) cr + mrj 1 o 2 , 

from the conditions of the constraint conservation 

= {$ 2) Wi} = 0, (2.15) 

wher e { } are the generalized Poisson brackets, see 
Ref. (Hf. 

The presence of the second-class constraints in a sys- 
tem implies that in the equations of motion and in the 
quantization procedure, the generalized Poisson brack- 
ets must be replaced by the generalized Dirac brackets, 
defined as 

{F, G} D = {F, G} - {F, $ 4 } C V] {$„ G} , (2.16) 

where the matrix (Cy) has the components Cy = 
In our case, 



' ( • jW-y)- 



(2.17) 
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The nonzero Dirac brackets of the basis fields have the 
form 



{77(x,t),77*(y ; i)} D = -i<5 3 (x-y), 
Mx,f),7r(y,t)} D = <5 3 (x-y). 

Hamiltonian equations of motion have the form 

V = {v,K} D = ( crV ) ? ? - ma 2 rf, 

if = W,n} D = (<r T V) 77* + ma 2 r], 

n = {n, %} D = i (cr T V) rf + ima 2 r], 

7T*={7r*,W} D =0. 



(2.18) 



(2.19) 



Of course, one can see from Eq. (|2.19[) that the evolution 
equation for r\ coincides with Eq. (|2.4I) . 

Using Eq. (|2.12l) . it is convenient to rewrite the wave 
equations for 77 and 7r as follows: 



f) = (<rV)?7 + i<72TO7T, 
7T = (<T T V) 7T + ia 2 mrj, 



(2.20) 



where the operation of the complex conjugation is ex- 
cluded. 

The two-component Weyl spinor field 77 is fermionic. 
When doing the canonical quantization of such a field, 
one should replace 77 and 7r with the operators 77 — > 77 
and 7r — > 7r, and define the equal-time anticommutators 
for these operators via the corresponding Dirac brackets 
(see Ref. [l4|). The nonzero equal-time anticommutators 
for the basic Heisenberg operators are 

[fj(x,t),7t{y,t)] + =i {?7(x } t),7r(y,t)} D |^ = - i7r= . 

= W 3 (x-y), (2.21) 

where the fundamental Dirac bracket is given in 
Eq. (|2TT5|) . Equation (j2~T2l holds true for the corre- 
sponding operators, such that we can rewrite Eq. ([2~2lj) 
as follows: 



[r){t,x),T)Ht,yj] = <5 3 (x-y). 



(2.22) 



Constructing the Heisenberg operators fj(x) that 
satisfy Eq. (|2.4p and the anticommutation relation 
[Eq. (H22)], we obtain 



fj(x) = 



d 3 p 

(2tt) 3 / 2 



2£ 
rn 



a_W- 



E+\V 

rn 



-ipx 



E 



(2.23) 



where E = \/|p| : 



e -W2 cos g/2 
e W2 S i n 61/2 



i<£/2 



sin 0/2 
e 1 ^ 2 cos 6»/2 



are chiral amplitudes; the angles and 9 fix 
the direction of the particle momentum, p = 
|p| (cos 4> sin 9, sine/) sin 9, cos 0); and a±(p) and a±(p) are 
the creation and annihilation operators, 



K(k),a CT ,(k')] 



=<W<5 3 (k-k') 
=0, 

=0, a = ±. 



(2.24) 



Note that Eqs. (|2.23[) and (12.24[) coincide with analogous 
expressions derived in Ref. [12| using heuristic arguments. 

Using Eqs. (l2"T4l . (J2H3I), and (|2~24]) . we get the total 
energy of the massive Weyl field via the creation and 
annihilation operators as follows: 

E tot = J d 3 rH = J d 3 pE (ala_ + 

+ divergent terms, (2.25) 

where "divergent terms" contains <5 3 (0) and can be re- 
moved by the normal ordering of operators. This result 
confirms the chosen interpretation of the operators a± (p) 
and a±(p). 

Let us introduce two propagators of the massive Weyl 
field in vacuum as follows 



Siix-y) =i(0\T{ V (x)rf(y)}\0), 
S 2 (x-y)=-(0\T{r,(x)r 1 T (y)}\0). 



(2.26) 



Here |0) is the vacuum vector for the annihilation opera- 
tors a±(p), and T is the sign of the chronological order- 
ing. With the help of Eqs. (f2T23|) and (f2T24|) we can cast 
the functions Si. 2 into the following forms fl3 |: 



Si{x) 
S 2 (x) 



d 4 P 



-\px 



(2tt) 4 to 2 - p 2 
d 4 p 



if 



cr 2 m 



(2tt) 4 to 2 -p 2 



if 



(2.27) 



where e is a positive, infinitesimally small constant, and 
the matrices a M are defined in Eq. pi)) . 

For our purposes, it is convenient to relate the intro- 
duced propagators [Eq. (|2.27|l ] to the propagator of the 
massive Dirac field. To proceed with the problem, let us 
introduce the Hamiltonian Weyl field S that is the pair 
7r and 77, 



(2.28) 



Using Eq. (|2.20p , one can find the the Hamiltonian Weyl 
field 5 satisfies the equation KE — 0, where the 4x4 
matrix operator K reads 



K _( 9 t - (er T V) -ia 2 m 
-\a 2 m d t - (ffV) 



(2.29) 
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The causal propagator for the Hamiltonian field S is 

S c (x-y)=(0\T{E(x)Ei(y)}\0), 
KS c (x) =-6 i {x). (2.30) 

Introducing the function S c instead of S c , 

S c = -6 2 S*c6 27 °7 5 , 6 2 = diag(a 2 ,a ), 

we derive for this function the following equation: 

(ir^-mr 5 )£ c =6 4 (x), (2.31) 

where = r 5 7^ and T 5 = — i7 5 . The 7 matrices T n — 
(r^,r 5 ) satisfy the usual 7- matrix algebra 

[r fc , T n ] =2r/ fc ", Vkn = diag(l, -1, . . . , -1), 

fc,ra=0,...,3,5, (2.32) 

so that S c can be interpreted as the causal propagator of 
the massive Dirac field. 

The propagators of the massive Weyl field [Eq. (|2.27l) ] 
can be related to the causal propagator S c of the Hamil- 
tonian field H, and then to the causal propagator S c of 
the massive Dirac field. First, we write 



PlS c Pl = - -((to - 0-3) ® 5*1, 

PlS c Pr =^(a 1 - ia 2 ) <8 5 2 , (2.33) 
i(<7 - a 3 ) ® Sx =P L S c P R r , 

i( CTl - ia 2 ) ® 5 2 -F L 5 c P L r°6 2 , (2.34) 

where the chiral projectors Pl,r were defined in Eq. I|2.7|l . 
Then, one can derive from Eqs. (|2.33l) and (I2.34[) that 



5i = - ia^ 1 (V» - A, 



o-a, (2.35) 



where S, P, V^ t , A M , and are the scalar, pseudoscalar, 
vector, axial-vector and tensor coefficients in the expan- 
sion of S c in the independent 7-matrix basis, 

S c = Sl + P7 5 + V M -f + A^V + T^ v , (2.36) 

where I is the unit 4x4 matrix and = (i/2)[7 M , 7"]-. 
Finally, using Eq. (|2.35l) we can express S1.2 as follows: 



Si 



-tr 



S-2 = — tr 



7 m7 °(1 - 7 5 )£ 



(2.37) 
(2.38) 



In Sec. IIV| we shall give the expression for S c in terms of 
the path integral. 



III. INTERACTION OF FERMIONS WITH 
EXTERNAL FIELDS 

In this section we briefly review the interaction of 
fermionic fields with external fields. We discuss the gen- 
eral situation when an arbitrary number of spinor fields 
are present. The cases of Dirac and Majorana fermions 
are compared. We also consider the gravitational inter- 
action of spinor particles. 

The most general classical Lagrangian describing the 
interaction of 1/2-spin fermions ip a , a = 1,2, ... , which 
are Grassmann- valued spinors, with a set of external 
fields has the form 15], 

-Lint ^Sab^a^b + ^ab^al^^Pb 

+ Vab^^b + A^a-f^lpb 

+ TjTab^aV^lpb + ^ab^aCT^tpb, (3.1) 

where s ab , 7r afc , V£ b , A^ fc , T£ b , and IP^ are the scalar, 
pseudoscalar, vector, axial- vector, tensor, and pseudoten- 
sor fields, respectively. Note that these external fields can 
depend on spatial coordinates. 

If we deal with Dirac fermions, the spinors %p a and 
■ipa are independent degrees of freedom. In this case the 
external fields in Eq. (13.11) arc Hermitian matrices in the 
indices a and b. 

If we have Majorana fermions, the spinors ?p a and ij) a 
are no longer independent. Instead they obey the rela- 
tion, cf. Sec, nn 



-0Q = Clp a , 



4>a = tC, 



(3.2) 



where C is the charge conjugation matrix which has the 
following properties (TTj : 



c = -c- 1 = -C T = -C\ 
c^c- 1 = -7^ T , c 7 5 cr 1 = 7 5T , 



C7VC" 1 = (7V) 



5\T 



n n-\ _ _ T 



(3.3) 



which are independent on 7-matrix representations. 

Using Eqs. (|3.ip - (|3.3l) . we immediately find that for 
Majorana fermions the matrices s a b, 7r a h, A^ b should be 
symmetric, whereas V^ b , T^ b , and U^ b are antisymmet- 
ric. This fact, in particular, means that the diagonal 
(i.e., when a — b) electromagnetic interaction of Majo- 
rana neutrinos vanishes (see also Sec. HVI below). 

The interaction of a spinor particle with a gravitational 
field can be implemented in a locally Minkowski space 
by the introduction of the covariant derivative there (see, 
e.g., Ref. |l6(). This covariant derivative also includes the 
contribution of the torsion, which, in principle, can have 
a nonzero value. The Dirac equation in flat space-time 
in the presence of torsion was also studied in Ref. (l7j . 
It was shown there that the contribution of the torsion 
is equivalent to that ofan external axial-vector field [see 
Eq. (13. ip ]. The issue of whether there is a difference 
between the gravitational interaction of Dirac and Majo- 
rana fermions was discussed in Ref. 181. 
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IV. MASSIVE WEYL NEUTRINOS IN 
EXTERNAL FIELDS 

Here we generalize the results of the previous section 
to include the interaction of the massive Weyl particles 
with external backgrounds. We discuss a particular case 
of massive neutrinos interacting with background matter 
and an external electromagnetic field. 

Note that in Sec. |TT] we discussed the case of a single 
massive Weyl field. It was also mentioned that a gener- 
alization to several mass eigenstates rj a with masses m a 
is not so difficult (see, e.g., Sec. IIII[) . The most general 
Lagrangian which involves the interaction of the fields 
r\ a with background matter and the electromagnetic field 
(E,B) has the form (see, e.g., Refs. 0,0), 

£ =Wa^d^ri a - ^m a (r]J<r 2 r] a - vl°2V*a) 



1 



+ O/t^J^B + iE)^] 
-^[eob^o-(E + iB)oi2»jS 
+ (£ t W?7> 2 (E-iB) C r7 ?b ], 



(4.1) 



where {[i a b) and (s a b) are the matrices of magnetic and 
electric dipole moments. 

The interaction with background matter is character- 
ized by the quantities g^ b , which are Lorentz four- vectors 
and Hermitian matrices in the particle species space, 



9 at — 9ba- A s was shown in Ref. [21j (see also Sec. IIII|) . 
the Hermitian matrices (fi a b) and (s a b) must be antisym- 
metric. 

The physical implementation of the Lagrangian in 
Eq. (|4.ip is a system of massive Majorana neutrinos prop- 
agating in background matter under the influence of an 
external electromagnetic field. In case we study the sys- 
tem of massive neutrinos, the components of the ma- 
trix (g^ b ) are related to the effective potentials of the 
active flavor neutrinos' interaction with background mat- 
ter, U» = (t/°,U) = diagfe,J7^,t/£), by means of 

the matrix transformation, (g^ b ) — U^U^U, where U is 
the mixing matrix. The zeroth component, U , is pro- 
portional to the effective number density of background 
fcrmions, and the spatial components, U, are the linear 
combination of the matter velocity and polarization. The 
explicit form of If 1 can be found in Ref. (l5j . 

The nonperturbative analysis of the complete sys- 
tem (|4.1[) which includes the nonzero, nondiagonal inter- 
action with the background matter and electromagnetic 
field, is difficult to perform. Nevertheless, it can be car- 
ried out in frames of the perturbation theory. Thus, first, 
we have to study the dynamics which results from the di- 
agonal terms of the Lagrangian in Eq. (|4.1I) . For the sake 
of the simplification of notations, we shall omit the index 
a: rj = rj a , = g% a , etc. Using the results of Sec. [Til wc 



obtain for the new fields an analog of wave equations in 
Eq. d23CT : 



v^ri - imcr 2 7r + ig^a^r/ = 0, 
cr*<9 M 7r — ima-2'n ~ ^9^ a ^i n = 0- 



(4.2) 



Similarly to Sec. [Til we introduce the causal propaga- 
tors S c {x-y) = (0|T{S(x)S%)}|0), where S T = (tt,T)), 
and the modified propagator S c = — &2S c &2j "f 5 ■ Using 
Eq. (|4.2I) . we get the equation for the modified propaga- 
tor S c in the following form: 



[n 1 (i^ + i 5M r 5 )-mr 5 ]5 c = <5 4 ( a 



(4.3) 



Note that Eq. (|4.3p is analogous to one obtained in 
Ref. [It} for the propagator of a massive Dirac particle 
interacting with an effective torsion field S£ s if we make 
the replacement g^ — > —S^. 

Then the path integral representation for S c can be 
found using the results of Ref. Such a representation 
has the form 



exp iT" 



d6 n 



(4.4) 



9=0 



where 9 n are Grassmann-odd variables that anticommute 
with T™, and 



Z = 



de / dxo / M(e)Be 

J Jen 



' <'a 

Bx I 



Bv \ Dip 

x exp {i (S* c i + Sgf) + ^»(1>"(0)} • 



(4.5) 



Here 



S c 



z 

Ye 



M 2 + x u d" 



+ ix ( + -^*d A 



dr, 



M 2 =m 2 + g 2 + 16^.9^ VW, 



(4.6) 



where e a p\a is the totally antisymmetric tensor with 
£0123 = 1) the values ip n are Grassmann-odd variables, 
the term 



<Sgf 



[fre + vy\ dr 



(4.7) 



is the gauge-fixing action, and the measure Ai{e) reads 



M(e) 



£>pexp 



ep dr 



(4.8) 



In Eqs. (|4.5I) - (|4.8[) . x(r), p(r), and 7r(r) are even trajec- 
tories, whereas x( T ) and v{t) are odd ones and e is the 
even variable. 
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One can interpret S c \ as the parametrization-invariant 
pseudoclassical action of a massive neutrino moving in 
the background matter. It should be noted that the first 
example of such action was the pseudoclassical action of 
the massive Dirac particle in 3 + 1 dimensions presented 
by Berezin and Marinov in Rcf. [22]. 

It should be noted that the path integral representation 
of the propagator in Eqs. (|4.4[) - (I4.8[) could be very useful 
in calculating the propagator of a massive Weyl field in 
background matter since it accounts for all the loop cor- 
rections. Moreover, here we do not restrict ourselves to 
the studies of homogeneous matter with characteristics 
like number density constant in time, as was assumed in 
Refs. dHH]. 

It is interesting to discuss the limiting case m — > 
of a massless Weyl particle propagating in background 
matter. Considering such a limit in Eq. (|4.6[) . we obtain 
the following pseudoclassical action: 



^cl — 



2 e (^A 1 ec ^ 



i 



dr. 



(4.9) 



Let us compare Eq. (14.91) with the pseudoclassical ac- 
tion of a massless Weyl particle proposed in Ref. 25;] : 



s (o) 



dr, 



(4.10) 



where a and 6 M are even variables. The quantization of 
the action in Eq. (|4.10p reproduces the quantum theory 
of a Weyl particle. Namely, one gets the description of 
a right (a = 1) or left (a = —1) massless neutrino in 
terms of a bispinor ^q(x), which satisfies both the Dirac 
equation with zero mass and the Weyl condition 



id^ (x) = 0, 



( 7 5 -a)* (x) =0. (4.11) 



On the contrary, the quantization of the action in 
Eq. (|4.9[) would recover the description of a massless 1/2- 
spin fermion in terms of two-component Weyl spinors t]q 
[see Eqs. (|2.9I) and (|2.10p ]. Note that some alternative 
pseudoclassical actions for massless fermions were dis- 
cussed in Rcf. 26]. 

Let us choose in Eq. (14. 10)) the gauge where a = — 1 
(left neutrinos) and 6 M = 2ieg Pl . Accounting for the fact 
that ip 5 =const if m = 0, we can cast Eq. (|4.9[) in the 



form S, 



(0) 



S- 



(0) 



(0) 



Sj° t , where 



- 8e W 



(4.12) 



can be regarded as the pseudoclassical action that de- 
scribes the interaction of massless Weyl neutrinos with 
the background matter. 

Finally, let us discuss how to find path integral repre- 
sentations for the propagators of a massive Weyl field 
5i.2, defined in Sec. |TTJ For this purpose we can use 
Eq. (|2.35[) . Let us decompose the generating function 
Z in Eq. (j4T2"|) in the variables 9: 



Z = J2 



Ho— in 



(4.13) 



ri=0 



where fi ...i n are the coefficients, which arc antisymmet- 
ric tensors at n > 1. For n = 0, /o = const, and n = 1, 
the quantity f n is a vector in a five-dimensional space. 

Using Eqs. (|4.4I) and f|4. 13[) . we represent S c in the 
following form: 

71=0 



1 



(4.14) 



On the basis of Eq. (|4.14l) , we can find the coefficients in 
the expansion of S c [Eq. (|2 .36[) ] as follows: 



S —fa — i/oi235, P — /5 + -^ £flIyX(J fftuXai 



V — - f r ^^afSXa r 



Au —fu 



T af3\cr 



<Jfi i 



Tfiv — 2 \Jliv5 J j 



(4.15) 



where we have used the representations of the propaga- 
tors Si t 2 in Eq. (|2.35[) . Here r\^ LV = diag(l, —1, —1, —1) is 
the metric tensor in the Minkowski space. Then, on the 
basis of Eqs. (|2.35p and (14. 15)) we can get path integral 
representations for S\^.- 



V. MASSIVE DIRAC NEUTRINOS IN 
EXTERNAL FIELDS 

In Sees. Inland IIV( we studied both free massive Weyl 
field and the interaction of such a field with different 
backgrounds. The most probable implementations of 
such systems are massive Major ana neutrinos. However, 
as was mentioned in SecU one cannot exclude the possi- 
bility that neutrinos are Dirac particles. That is why we 
consider below neutral massive Dirac particles possess- 
ing anomalous magnetic moments, which interact with 
background matter and with an external electromagnetic 
field. 
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Using the results of Sec. IIII1 we can write down the 
most general Lagrangian which describes the interaction 
of several massive Dirac particles ip a with matter and an 
electromagnetic field. On the basis of Eq. p. II) we should 
identify T^ b with n a j ) F^ 1 ', where (fJ, a b) are anomalous 
magnetic moments and F^ v = (E, B) is the tensor of the 
electromagnetic field. The interaction with background 
matter can be introduced if we set V£ b = — A^ b = g^ b /2 
in Eq. (|3.1I) . where the quantities g£ b have were already 
introduced in Eq. (|4.1|) . Finally, we get this Lagrangian 
in the following form (see also Ref. |27| and references 
therein) : 



(5-1) 



where ^ = 7^(1 - 7 5 )/2 . 

As we already mentioned in Sec. lIIIl the matrix of mag- 
netic moments, (fi a b), is Hermitian for the Dirac neutri- 
nos. It means that nonzero, diagonal elements of this 
matrix are possible 28]. The diagonal element /i aa cor- 
responds to a magnetic moment of the mass eigenstate 

The exact analysis of the dynamics of the system de- 
scribed by the Lagrangian in Eq. (|5.1[) with nonzero, non- 
diagonal elements of the matrices {g^ b ) and (/x a &) is com- 
plicated but can be carried out in frames of the pertur- 
bation theory. That is why, as in Sec. IIV1 we will discuss 
only diagonal terms in Eq. (|5.1[) with a — b. In this case, 
we get the wave equation for the propagator of a single 
Dirac neutrino S^: 



1 



i 7 ^-m--y(l- 7 5 )«? 



-a F^ 



Sn(x) = S 4 (x). (5.2) 



where m = m a , g^ = g% a , and /i = fx aa . Then, using 
Eq. (|5.2p . we obtain the following equation for the mod- 
ified propagator, Su — Sbr 5 : 



-r 5 g„ 



mF 5 



I^t^y" f^ 



S D (x)=6 4 (x), (5.3) 



where the matrices F M and L 5 were defined in Sec. [TTJ 

One can see that the vector part of the interaction with 
matter is equivalent to the interaction with an effective 



electromagnetic field: (qA^ 



(7 M /2, where q is an ef- 



fective electric charge. The axial vector part of the inter- 
action with matter is equivalent to the interaction with 
an effective torsion field: 5^ ff = —g^/2. Contributions 
due to such interactions to the propagator of the Dirac 
massive particle were studied in Ref. [l7| . The analysis of 
the propagator of a fermion with an anomalous magnetic 
moment was made in Ref. 12911 . 



Using the results of the latter works, we can obtain the 
path integral representation for the propagator Sd . The 
general structure of such a representation is similar to 
one given by Eqs. P~4"|) . (|4T71) . and (O with the 

substitution of the pseudoclassical action in Eq. (I4.6P by 
the following one: 



S c 



+ \x ( K,^ + 



dr, 



(5.4) 



where 



Ml=ml 



d^g„ - d v g^, 



and mo = m — 2ifiF a pip a ip" . 

Again, we can mention that Eqs. (|4.4[) . (|4.5I) . and (|5.4[) 
are the most complete expression for a propagator of a 
massive neutral Dirac particle with an anomalous mag- 
netic moment, since they exactly account for the influ- 
ence of background matter and an external electromag- 
netic field. Besides the loop corrections, these expressions 
also contain the inhomogeneous matter contributions, 
since we do not suppose that d^g^ = 0. Thus, Eqs. (14.51) . 
(14~4|) . and ((5^4)) generalize the results of Refs. |23l. |24|. 
where the neutrino propagator in homogeneous matter 
was obtained on tree level without loop corrections. 



VI. CONCLUSION 

In the present work, we have studied massive 1/2- 
spin particles in various external backgrounds, keeping 
in mind applications to neutrino physics. We have been 
mainly interested in massive Majorana (or Weyl) fields. 
However, massive neutral Dirac particles have also been 
considered. We have formulated classical Lagrangian the- 
ory of the massive Weyl field in terms of Grassmann-odd 
two-component spinors. Then, we have constructed the 
Hamiltonian formulation of such a theory, which turns 
out to be a theory with second-class constraints. Us- 
ing this formulation, we have canonically quantized the 
massive free Weyl field. We have derived propagators of 
the Weyl field and related them to the propagator of a 
massive Dirac particle. 

Then we have studied the massive Weyl particles prop- 
agating in a background mater. We have found the path 
integral representation for the propagator of the massive 
Weyl field in background matter, as well as we have ob- 
tained the corresponding pseudoclassical action for mas- 
sive Weyl particles. The massless limit of such an action 
was compared with the results of other works. Finally, 
we have studied the path integral representations for a 
neutral massive Dirac particle with an anomalous mag- 
netic moment moving in background matter under the 
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influence of an external electromagnetic field. From this 
representation, we have derived the pseudoclassical ac- 
tion of the corresponding neutral massive Dirac neutrino 
with an anomalous magnetic moment. 

The results of the present work can be applied for the 
study of the propagation of massive mixed (Majorana or 
Dirac) neutrinos in dense matter and strong electromag- 
netic fields. As demonstrated in Ref. [30(, path integral 
representations of particle propagators allow one to effec- 
tively calculate the propagators. We hope that the repre- 
sentations derived in the present work will be helpful for 
describing the neutrino motion in various astrophysical 
and cosmological media (see, e.g., Ref. [HJ). 

In the case of a Dirac neutrino moving in a singular ex- 
ternal background one should carefully define boundary 
conditions for corresponding solutions, as this is equiv- 
alent to choosing a self-adjoint Hamiltonian in such a 
background. To this end, one can use a general the- 
ory of self-adjoint extensions of symmetric operators and 
some analogy with the motion of charged particles in an 
Aharonov-Bohm field (see, e.g., Ref. |32j|). 

It should be noted that neutrinos may have magnetic 
moments (for a review, see Refs. [13, Hll)- The best ex- 
perimental constraint on the Dirac magnetic moment of 
an electron neutrino is ^„ e < 2.9 x 10 _11 /zb [34 |. where 
/is = e/2m e is the Bohr magneton. Astrophysical con- 
straints on the Dirac magnetic moments are stronger: 
li v < 1.1 x 10 _12 /jb [35j . Note that in order to sat- 



isfy m u < 1 eV, there should be a more natural scale 
for the Dirac magnetic moments of a neutrino: \i v < 
10~ 14 /iB [36j ] . If neutrinos have magnetic moments, an in- 
homogeneous magnetic field, acting on such particles, can 
change their kinetic energy, i.e. produce work. This fact 
implies the possibility of the neutrino creation form the 
vacuum by strong inhomogeneous magnetic fields (see, 
e.g., Ref. [37j). In the case of Dirac neutrinos, this effects 
was studied in Ref. (38|. In this connection, it is inter- 
esting to generalize the technique of the present work to 
quantize Majorana neutrinos that have anomalous transi- 
tion magnetic moments, and then to study the creation of 
such neutrinos by strong inhomogeneous magnetic fields. 
In addition, by applying similar methods, one can study 
the creation of Majorana neutrinos from vacuum by the 
inhomogeneous background matter. 
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